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We simultaneously analyze vector and axial-vector spectral functions in vacuum using hadronic 
models constrained by experimental data and the requirement that Weinberg-type sum rules are 
satisfied. Upon explicit inclusion of an excited vector state, viz. p' , and the requirement that the 
perturbative continua are degenerate in vector and axial-vector channels, we deduce the existence 
of an excited axial- vector resonance state, a'l, in order that the Weinberg sum rules are satisfied. 
The resulting spectral functions are further tested with QCD sum rules. 

o 

(N 

. Chiral symmetry is believed to be a fundamental symmetry of QCD. However, this symmetry is spontaneously 
broken at low temperatures and chemical potentials by a finite expectation value of the quark (or chiral) condensate, 
, (0|gg|0) ~ —2 fm~^ per light quark flavor in vacuum. At higher temperatures, the symmetry is expected to be restored 
as the value of the condensate approaches zero. Such behavior has been observed in lattice-QCD computations [1, 2], 
] but it remains a long-standing goal to observe chiral restoration in experiment. This is difficult because a direct 
C^' measurement of the chiral condensate is not possible. 

^' One way to indirectly infer the chiral condensate and deduce chiral symmetry restoration is through the use of 
I sum rules. They are formulated to relate hadronic spectral functions to properties of the ground state (condensates) 
' including chiral order parameters. The observation of in-medium changes of hadronic spectral functions can then 
signal pertinent changes in these order parameters. Of particular interest are spectral functions of so-called chiral 
' partners, i.e., hadronic states which are degenerate at chiral restoration, but are split in vacuum due to dynamical 
^ , chiral breaking. The prime example are the iso- vector vector and axial-vector channels which are connected through 
Q^ ■ Weinberg-type sum rules 0-01 (see, e.g,, Ref. for a recent study of chiral partners in the open-charm sector). 
' Therefore an accurate measurement of both the vector and axial- vector spectral functions can be used to infer chiral 
I order parameters and thus "observe" chiral symmetry restoration. Additional information can be gleaned from the 
■ spectral functions by considering QCD sum rules 0, which are dependent on both chirally symmetric and chirally 
, breaking ground-state properties. A lot has been learned about the in-medium vector spectral function from dilepton 
' measurements and their interpretation but no experimental measurement of the in-medium axial- vector spectral 

function has been performed to date. One is thus left with constructing effective models of the axial- vector spectral 
' . l function in order to study chiral symmetry restoration. 

I> In order to reliably evaluate in-medium effects on the vector and axial-vector spectral functions, one should first 

control their properties in vacuum. This is aided by accurate experimental measurements of the vacuum spectral 
functions through r-decays by the ALEPH and OPAL Collaboration s flll[T2j . On the one hand, sum rules have been 

d . used along with experimental data to calculate the condensates |13l - [la |. On the other hand, the condensates can be 
used to constrain the spectral functions. Weinberg's original study assumed a J-function approximation for both the 
low-lying vector and axial-vector resonances in his sum rules to determine the famous relation between the masses 
of the p and ai mesons, nia^ = V2mp [3. Similarly, QCD sum rules have been applied in the vector channel, where 
early works approximated the p by a (5 function while also considering medium effects [l6l . [iTj . More recent work has 
considered spectral functions with a Breit-Wigner shape for the p peak [l^ , and have included the axial- vector channel 
in their analyses [l9l-l2l|. while others analyzed calculations based upon a microscopic theory [20l - [23j . Clearly, the 
focus in previous work was on QCD sum rules, with the contribution to the continuum approximated by a 0-function 
with a threshold energy. Furthermore, until the most recent studies [13, [2l|, the vacuum spectral functions did not 
take advantage of the high-precision r-decay data. The often simplified constructions of spectral functions used to 
analyze the sum rules, and the dearth of analyses to include the axial-vector channel or consider the Weinberg-type 
sum rules, leaves this area open for further considerations. 
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In the present paper, we simultaneously analyze vector and axial-vector spectral functions in vacuum using an 
extended model which combines a microscopic p spectral function with Breit-Wigner ansatze for the ai and the first 
excited states. This model is quantitatively constrained by both the experimental r-decays [HI, and the Weinberg- 
type sum rules [B^Hl • By using this combination of criteria, the model may be considered a non-trivial fit of the data. 
Novel features of our analysis include the study of excited states and the postulate that the continuum contribution is 
identical for both the vector and axial- vector channels, as should be the case in the perturbative regime. In particular, 
the use of the Weinberg-type sum rules leads us to deduce the presence of an excited axial-vector state, about which 
rather little is known to date 24] . As an additional check, we utilize the constructed spectral functions in a pertinent 
analysis of QCD sum rules. 

The outline of this paper is as follows. In Sec. [Ill the sum rules used in our investigation are introduced. In Sec. lIIIl 
we detail the main ingredients to the vector and axial- vector spectral functions. Section ITVl presents the results of our 
fitting procedure and discusses our main findings in the context of Weinberg sum rules. In Sec. |Vl the constructed 
spectral functions are implemented into QCD sum rules, and we conclude in Sec. IVIl 



II. WEINBERG AND QCD SUM RULES 

Sum rules are a valuable tool for understanding non-perturbative aspects of QCD as they relate two different 
formulations of the same correlation function to each other. For the current paper, we will focus on two classes of 
sum rules, Weinberg-type and QCD sum rules. 

Let us begin by defining the current-current correlators for the vector and axial- vector channels, 

nr(9) = d^xe^^^Tfy{x)fym, (1) 

nr(9) = d^xe^^-(r4(x)jl(0)), (2) 

where jy = i {u'y^u — dj^d^ and — ^ (uj^j^u — d'j'^j^d^ are the pertinent currents in the quark basis. The 
correlators can be decomposed into 4-dimensional transverse and longitudinal parts as 

K':Aiq') - [-9'^ + ^) n^.A('z^) + '-^^U^') ■ (3) 

Since the vector current is conserved, Ily = 0. The longitudinal component of the axial-vector is governed by the 
contribution from the pion which causes the axial current not to be conserved. In vacuum, the self-energy of the pion 
is assumed to be negligible, rendering the imaginary part of the longitudinal axial- vector polarization of the simple 
form 

Imni(s) - -7:f,sS {s - ml) (4) 

It is also useful to define an additional polarization function, 11^(9^) = nj(q'^) +ll^{q'^) (which corresponds to 11^ 
in Ref. [^). The transverse polarization functions for both vector and axial- vector channels are used to define the 
spectral functions such that 

pv{q') ^-ilmn^(g2), (5) 
PAiq') ^ -ilmn5(q2) . (6) 

A spectral function for can be defined in a similar manner as 

PAiq^) = --ImlLAiq^). (7) 

TT 

From the definition of fl^ and Eqs. dH) and (O, we see that pa{s) = pa{s) + f^s5 (s — mj), i.e., this spectral function 
has contributions from both the pion and the axial- vector mesons (plus continuum, see below). 

The Weinberg-type sum rules characterize moments of the difference between the vector and axial-vector spectral 
functions, thereby quantifying the effects of chiral symmetry breaking. They have the general form 



/"OO 

/ dss"Ap(s) = /„ , 

"'0 



(8) 
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where Ap = pv — Pa-, n is an integer, and /„ are chiral order parameters dependent on the sum rule in question. 

The first two of these sum rules were originally derived by Weinberg Q (hence the name for the class of sum rules) 
using current algebra arguments. They read 

(WSRl) rds^ = fl, (9) 
Jo s 

(WSR2) / dsApis) = f^ml = -2mg{qq) , (10) 
Jo 

and correspond to n = —1 and n = in Eq. ([5]). Here, ~ 5 MeV refers to the average current light-quark mass, 
while the Gellmann-Oaks-Renner relation [25], 

/>2 ^ -2m,{qq), (11) 

was used to obtain the second equality in Eq. (|10p . In Weinberg's original work, the chiral limit was considered; 
here we have included terms linear in the quark mass 26-29j. Corrections to higher powers in uiq [26-28, 30] are 
expected to be small. ^ One sees that the chiral condensate can be calculated using Eq. (fTO)) should both the vector 
and axial-vector spectral functions be known precisely. The sum rule for n = —2 was introduced by Das, Mathur, 
and Okubo H, 



(WSRO) / ds^^-f^rD-FA, (12) 



s2 3 

where (rj) is the mean squared radius of the charged pion, and Fa is the coupling constant for the radiative pion 
decay, ■n'^ — )• /i^t'p7. We label this sum rule as the 0**^ one since it has a smaller value of n as compared to Weinberg's 
initial sum rules. Lastly, Kapusta and Shuryak [sj derived a sum rule for n — 1, 

/•oo 

(WSR3) / dssAp(s) = -27ra,(C'4) , (13) 

where (©4) is the part of the four-quark condensate which breaks chiral symmetry. The explicit quark content of this 
operator is given by 

(O4) - \ (((?I7m75A"?/ - d7^75A'^d)' - {u^^^X'u - d7^A'^d)')) , (14) 

where A° denote the Gell-Mann matrices. It is common to assume that this operator can be factorized into the chiral 
condensate, 

1 fi 

(O4) = -n{qqf , (15) 

where k is a parameter larger than one to mimic the contributions beyond ground-state saturation. Typically the 
Wcinberg-type sum rules are expressed in vacuum in terms of pA as is done here. However, it is possible to include 
the pion pole into the axial-vector spectral function and thereby write the sum rules in terms of pA with their 
right- hand- side (RHS) appropriately adjusted. This was done in Ref. Q for the sum rules at finite temperature. 

For studies of chiral symmetry restoration, satisiying these sum rules is critical both in medium and in vacuum. 
In the present work, we will make the first step by using the first three, Eqs. (O, pUj) . and p^ . to constrain the 
parameters of the vacuum spectral functions. The last sum rule, Eq. (I13p . is not used due to a large uncertainty in 
the value of k and the rather high sensitivity of the integration to large s where the control from experimental data 
is limited. Nevertheless, we will still examine how well it is satisfied by the constructed spectral functions. 

The QCD sum rules, on the other hand, apply to the vector and axial-vector channels separately. They were first 
introduced by Shifman, Vainshtein and Zakharov using the operator product expansion (OPE) 0, Hj . They relate the 
integral over the spectral function to a series of ground-state operators (expectation values) specific for the channel 



^ Note that the sum rule in Ref. [s^ labeled as Weinberg's 2"^* sum rule is not the same as the one considered here; the one considered 
here corresponds to the convergent sum rule in the equal quark-mass limit of Ref. [30| . 
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of interest. To improve the convergence of the integral over the spectral function, one perforins a Borel transform on 
both sides of the sum rule. For the vector channel one obtains 

M2 S ^ " 87r2 + TT J + Af4 + 24M4 ^ TT ^^-^ 81M6 ) ■ • • , (if J 

and for the axial-vector channel 

^ rds^e-^/''^ = J_ fl + "^^^^"'^^ + -^(^G' ) + ) (17) 

ApJo s 87r2 + TT J + M4 +24Af4^^^M-^+ 8i^^6^^4;--- • U^J 

By performing the Borel transform, one trades the space-like 4-momentum, — —Q^, with the Borel mass, M^. 
Note that the axial- vector spectral function is defined to contain the contribution from the pion pole (use of pA, not 
Pa)- To the order which we are working, the following operators figure: the chiral condensate, (qq), the gluon 

condensate, (^G^^), and the vector and axial-vector four-quark condensates, (OY) and (O^), respectively. The 
four-quark condensates can be expressed in terms of their quark content as in [1^ 

(Ol) = ^ ((S7p75A'^" - d7p75A'^d)') + ( (^7mA"" + dj^.^'^d) 97^A'^0 - (18) 

\ q—u,d,s I 



\ q—u,d^s I 

and are related to the chirally breaking four-quark condensate of Eq. as 

16/7 ,,^y^ , 11 , 



(04) = y (^^(Ol) + ^(Of>j- (20) 

These forms of the four-quark condensates and the coefhcients in the sum rules are chosen such that (Ojf) and {O^) 
have a simple factorized form, viz. {OX) = Kv{qq}'^ and {O:^) = ka{ciq}^- The parameters kv and ka are, in principle, 
independent and thus could take on different values; they are related to k in Eq. (IT5|) via 

K — — Kv H Ka ■ (21 

18 18 ^ ^ 

For simplicity we have chosen Ky and ka as numerically identical, which then implies from Eq. (|2ip that k is also 
numerically the same, viz. kv — ka — k. A more general discussion on the properties of four-quark condensates can 
be found in ^sT]. We note that the sign of the chiral-condensate term is the same for the two sum rules, not opposite 
as one might expect since the chiral condensate is a chirally odd operator. The latter reasoning actually applies to 
the transverse part of the axial-vector current. However, the QCD sum rule written above is for the total current 
including the longitudinal part induced by the pion contribution. This leads to the signs presented in Eq. (jl7p . We 
also note that the Weinberg sum rules 1-3 can be derived from the QCD sum rules by subtracting both sides of the 
vector and axial-vector QCD sum rules from each other, Taylor expanding the Borel convergence factor e~*/*^ , and 
equating the coefficients of equal powers of on each side. A similar procedure was pointed out in Ref. @ and used 
in Ref. [5] to derive Eq. p^ . 

If one knows the vector or axial-vector spectral function, the QCD sum rule can be used to determine the values 
of the condensates. Conversely, since the condensates are universal low-energy operators, they can be calculated 
from other processes, and in turn be used to constrain the spectral functions. Though related to the Weinberg-type, 
QCD sum rules can provide different constraints on the spectral functions because both chirally even and chirally 
odd operators are involved. A recent study of the constraints on the spectral functions imposed by the chiral odd 
operators in the QCD sum rules can be found in Ref. [s^- A potential drawback of the QCD sum rules comes from 
the relatively large uncertainty of the values of the gluon and 4-quark condensates, as compared with the pion mass 
or its decay constant in the first three Weinberg sum rules. Therefore, we will not use the QCD sum rules to constrain 
the vacuum spectral functions, but rather use the constructed spectral functions to constrain the condensates and in 
this way perform a consistency check. 

In principle, there are corrections to the QCD sum rules 0, H, H^- Besides the conventional perturbative and power 
corrections, the most prominent one not included here is associated with the exchange of instantons with the vacuum. 



It turns out that the importance of instantons depends on the specific correlator [33[, i.e., whether direct-instanton 
interactions are operative, like in the scalar and pseudoscalar meson channels. In both the vector and axial-vector 
channels these are absent and the remaining corrections are probably small [11] . We therefore believe that the QCD 
sum rules considered here are sufhciently accurate. 
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III. HADRONIC MODELS FOR VECTOR AND AXIAL- VECTOR SPECTRAL FUNCTIONS 



So far, we have discussed Weinberg-type and QCD sum rules and their usefuhiess in constructing vacuum spectral 
functions. In this section, we present a model which we will use to construct the spectral functions suitable for a 
physically motivated fit to r-decay data. Note that for the purpose of the present work the concrete fit functions and 
parameter values are not important as long as the data are accurately reproduced. 

Our main ansatz, which is one of the differences of our work from previous analyses, is that the vector and axial- 
vector spectral functions are divided up into three parts: the ground-state resonance, a first excited state, and a 
universal continuum, 

pvis) - p^^is) + p'^is) + (22) 
Pa{s) = pf(s)-KpX(s)+P™"*(s)- (23) 

The explicit form of each part will be discussed in turn. 

The ground-state resonance in the vector channel, the p(770), has been well-studied in effective hadronic La- 
grangians; we here employ the spectral function of Ref. 34] which was originally fit to the pion electromagnetic 
form-factor and P-wave tttt scattering phase shifts, but also turns out to describe the experimental r-decay data 
well [ssj; its medium modifications have been widely ap plie d to experiment, e.g., in dilepton and photon production 
in heavy- ion collisions and in elementary reactions [37[. These features make the p spectral function a suitable 
starting point (both in vacuum and for future studies in medium) for the other two components which will largely 
rely on fits to the r-decay data. 

The spectral properties of the ai meson are much less studied, especially in the medium. Since our eventual 
objective are studies of the spectral function's medium modifications (and the pattern of chiral symmetry restoration 
for the vector and axial- vector), we adopt a more schematic ansatz for the ai spectral function. Rather than using the 
vacuum spectral function from an effective hadronic model (see, e.g., Refs. |38l - l44 |). we employ a generic Breit-Wigner 
form, 

Pa, (S) — ^ n , (24) 

where Ma^ is the bare mass, gai the axial- vector coupling constant, and (s) the energy dependent width of the 
resonance. These parameters will be determined uniquely for the oi. The form of the coefficient M^_^/g'^_^ is the same 
as in vector-meson dominance. The energy dependence of the ai width will be approximated by an 5*- wave decay 
into pn. In order to properly incorporate the three-pion final state, the spectral shape of the p will be accounted 
for by integrating over its spectral function. Finally, to simulate the finite size of the pnai vertex and to control the 
large energy behavior of the ai spectral function, a hadronic form-factor is included. Putting all of this together, the 
energy dependent width of the oi takes the form 



where -Dp(s) is the p propagator, and the center of mass momentum of the p and tt decay products, pj?^(s, s'), is given 
by 

PiUs, s')=^-\ ^" ' " . (26) 




This ansatz involves two additional parameters, a form-factor cut-off, Aq^ , and the pvrai coupling constant which 
determines the magnitude of the on-shell width, F^^, defined as F°^ — Fai(M^J. 

The first excited resonances in both vector and axial- vector channels, namely the p' and a'l, are included in the 
current construction. Their contribution to the spectral functions is presented here while the motivation will be 
provided later. Again, we note that not only is this the first attempt to include them in such a study, but given the 
well established p' contribution, we will be deducing the need to include the a'l based upon the Weinberg-type sum 
rules. A microscopic description of the excited resonance states is not well understood to date. Therefore, as was 
done for the oi, the spectral functions of the excited resonances will be approximated by a Breit-Wigner shape as in 
Eq. (|24p . To construct the energy-dependent width, we are once again guided by plausible decay channels. However, 
unlike for the ai case, the decay products are not well established. Instead we postulate that the width exhibits a 
threshold which is controlled by an estimated mass scale. We further postulate that the partial-wave distribution 
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of the decay products is the same as for its corresponding ground state, namely the p' decays through a P-wave 
process while the a\ decays through an S'-wave process. Furthermore, since the actual decay products of the excited 
resonances are not explicitly specified, treating them off-shell by folding in their spectral function is not warranted. 
As for the ground states, a form-factor will be included. This leads to a width of the p' as 



/ 



3/2 



th 



(27) 



and for the a'^ as 



1/2 



A? 



A' 



Ml 



(28) 



correspond to the 
Lastly, the overall 



The parameters Mpi and M^'^ set the bare mass for the pertinent resonances; ■* and M^^^' 
effective mass thresholds of the decay products, while Ap/ and A^'^ set the form-factor scales, 
strengths of the widths are given by F^, and F°, , which again have been normalized such that Fp/ ( Af^, ) = F^, and 
rai(Mf,) = FO,.^ 

At high energies, the contribution from the continuum can be calculated from perturbative QCD. Most previous 
work has assumed that this contribution exhibits a threshold energy below which it vanished while the perturbative 
value was assumed above. We have adopted a continuous function for all energies which approaches the perturbative 
value for large energies |5|], 



P 



1 



(29) 



The parameter Eth plays the role of a threshold energy while 5 determines how fast the limiting value is achieved. 
Furthermore, because at high energies, i.e., in the perturbative regime, QCD is chirally invariant, we postulate an 
identical continuum contribution to both channels. Thus the parameters Eth and 5 are the same for the vector and 
axial- vector channels and p™"* — p^f"'- This is the first time that this feature has been implemented; besides the 
underlying physical motivation of degenerate perturbative continua, it will play a crucial role in interpreting the 
vacuum spectral functions. 



IV. VACUUM SPECTRAL FUNCTIONS AND DISCUSSION 



The total number of fit parameters of our model for the vector and axial-vector spectral functions in vacuum is 
16 (the three from the p spectral function have been fixed before). They include the masses Mx, the couplings 
gx, the width strengths F^, and the form- factors scales Ax for each X — ai,p' and a[ states. The p' and a[ 

(X) 

furthermore incorporate a threshold mass scale, M^^ , and finally the continuum carries the two parameters Eth and 
S. The continuum also depends on the strong coupling constant which was chosen to be as(lGeV) = 0.5. This set 
of convenient and physically motivated parameters was determined by fitting the spectral distributions of hadronic 
T-decays into an even and odd number of pions (ill . Il2j . In addition, the requirement of reproducing the Weinberg- 
type sum rules was imposed (except the third one, Eq. (I13p . cf. our comment above). The fit was not done for all 
16 parameters simultaneously, but rather three sequential fits were performed. First, the parameters of the ai peak 
and the continuum were determined from the axial-vector ALEPH r-decay data. Second, the p' peak parameters 
were determined from the vector ALEPH r-decay data. And lastly, the parameters associated with the a'^ peak were 
determined from the Weinberg-type sum rules. The parameter values are summarized in Tab. HI Figure [1] shows the 
resulting spectral functions, divided by s to render them dimensionless, compared with the experimental data. A 
measure of the goodness of fit is the coefficient of determination, R^. The vector channel has an value of 0.991 
while the axial- vector channel (without the a[) has an R^ value of 0.997. Both of these indicate very good agreement 
of the resulting spectral functions and the experimental data. 

The results for the Weinberg-type sum rules are graphically assessed in Fig. [21 Here the left-hand-side (LHS) of 
each sum rule is plotted as a function of the upper limit of the energy integration. Toward high energies, i.e., when 
the spectral functions degenerate so that their difference no longer contributes to the sum rule, the curve should 
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ai 


P' 




Mx 


1.246 GeV 


1.565 GeV 


1.802 GeV 


^ X 


0.612 GeV 


0.32 GeV 


0.2 GeV 


ax 


6.15 


11.44 


28.70 


Ax 


0.61 GeV 


1.41 GeV 


1.24 GeV 




N/A (Stt) 


0.56 GeV 


0.96 GeV 



Continuum parameters 



Eth 
S 



1.60 GeV 
0.227 GeV 



Constant parameters 




139.6 MeV 


U 


92.4 MeV 




5 MeV 


Qs(lGeV) 


0.5 


(rl) 


0.439 fm^ 


Fa 


0.0058 


im) 


(-0.25GeV)3 



TABLE I: List of parameters of the constructed model as constrained by experimental data and the Weinberg-type sum rules. 




FIG. 1: Spectral functions for the vector (left) and axial-vector channel (right) compared to experimental data for hadronic r 
decays by the ALEPH collaboration [l3|. The different curves highlight the contributions to the total spectral function (solid 
curve) from the ground-state resonance (dotted curve), the excited resonance (dashed curve), and the continuum (dot-dashed 
curve) . 



converge to the value of the RHS, represented by the dashed curve. Table [TT] quantifies the numerical deviation of 
the asymptotic value from the RHS. The values are quoted such that a positive (negative) deviation means that the 
contribution from the vector (axial- vector) channel is too large. Primarily by introducing, and then adjusting the mass 
and the coupling of the a'l state, "perfect" numerical agreement with the RHS of Weinberg-type sum rules 1 and 2 can 
be achieved. This was intentionally done since the respective RHS of these two sum rules are known with much better 
precision than for WSR-0 and especially WSR-3. Not surprisingly then, the latter exhibits the largest deviation. The 
value of K which we use to quantify this deviation is actually determined from the QCD sum rules which will be 
described in the next section. In all cases, even for WSR-3, the deviation of the converged value from the asymptotic 
value dictated by the RHS is small compared to the size of the oscillations seen at low energies. The WSR-3 is more 
sensitive to the higher energy regime as compared to the other sum rules because of the larger power of s in the 
integral in Eq. (IT^ . Its deviation suggests that a little less spectral strength is needed in the axial- vector channel 
at high energies. Nonetheless, the "excess" axial-vector strength in WSR-3 is less than 10% of the total strength 
generated by the a^. Another possibility is the introduction of a second excited vector state at higher energies; we 
refrain from this possibility due to the lack of constraints available by the current scheme. 

A few comments pertaining to our fit are in order. First, because we have postulated the continuum for the 
vector and axial-vector spectral functions to be identical, the energy of the effective threshold is larger than in 
previous applications. This is essentially dictated by the dip region in the axial-vector spectral function data around 
s — 2.2 GeV^, which lies below the pQCD continuum level of (1 -I- ag/Tr) ~ 0.015. The onset of the continuum must 
therefore be pushed out to higher energies to accommodate this dip. This can be observed in Fig. 1(b) by inspecting 
the contribution from the continuum, the dot-dashed curve. The need for larger Eth may be somewhat mitigated by a 
tuning of 5. By allowing the onset to occur more slowly, i.e., with a larger value of (5, a slightly lower onset energy can 
be obtained. Nevertheless, even with this tuning, the resulting onset energy, Eth, remains significantly larger than in 
previous studies. Our value for the onset energy might rather be considered as a lower bound, since, in principle, the 
continuum could be pushed out to even higher energies and compensated by adding further spectral strength to the 
excited states (i.e., decrease gx or adding more states). However, a continuum threshold at even higher energies is 
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0.005 



-0.005 
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(c)WSR 2 
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(d)WSR 3 



So (GeV-) 



FIG. 2: Graphical representation of the Weinberg-type sum rules. The LHS of Eqs. (|9}, (|10p . (|12p and (|13p is plotted as 
a function of the upper integration limit (solid curve) , compared to the theoretical experimental value for the RHS (dashed 
curve). The dot-dashed curve is the same as the solid curve but excluding the contribution from the a'l state. 



WSR 






2nd 


3rd 


% agreement 


-1.28% 


-0% 


-0% 


-96% 



TABLE II: Percent disagreement between LHS and the RHS of the Weinberg-type sum rules resulting from our fit. 



not well constrained by data. Therefore we have taken the more conservative approach here. 

Second, the current study is the first to include excited resonance states into a construction of the spectral functions 
using sum rule techniques. Previous works [ToU^H have argued that their continuum ansatz, with a smaller threshold, 
was providing sufficient agreement with the data such that explicitly treating these states was not necessary. However, 
when the continuum is pushed to higher energies, a region in energy is created where the vector spectral function 
is no longer in good agreement with experiment (c/. Fig. 1(a) for the continuum contribution undershooting the 
data). To accommodate this, it is therefore natural to include the effects of a p' resonance. This argument, based 
upon observed bumps in experimental data, cannot be applied to the axial-vector channel, since there is no direct 
indication which suggests a clear need for an a'l state a priori. However, our enforcing of the Weinberg-type sum rules 
dictates additional strength in the axial-vector channel which is not incompatible with data. To further illustrate 
this point, let us switch off the a'l peak. The resulting axial- vector spectral function still describes the data as seen 
in Fig. [31 but evaluating the Weinberg-type sum rules in this scenario leads to appreciable discrepancies in all cases, 
cf. the dot-dashed curves in Fig. [51 Therefore, the requirement of satisfying the sum rules by the constructed spectral 
functions lets us deduce the presence of an a[ state. 

Third, one may ask if there is other evidence for an excited axial-vector state. The Particle Data Group 0] lists 
three possible candidates, ai(1640), ai(1930) and ai(2095), though none of them are very well established (the latter 
two can only be found in the "Further States" list and none are included in the summary table). Being the first 
excited resonance, it would be natural to associate our a[ with the ai(1640), but the mass and width do not match 
well. A similar problem actually holds for the excited vector resonance: while we fit the ALEPH data with a single p' 
state of mass 1565 MeV, the PDG has states at /ci(1450) and p(1700). This suggests that the bump seen in the data 
is an amalgamation of the two excited p states. In principle, we could have tried to use a spectral function with two 
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FIG. 3: Spectral function in the axial- vector channel without the a'l, compared to experimental data [Til ]. 

excited states, but the significance for the latter is not apparent from the single bump in the r-decay data and thus 
it would only serve to increase the number of our parameters without improving the fit. Though we call our state p' , 
it may well represent the contribution from the two excited resonance states quoted by the PDG. A similar argument 
holds for the a[ state. It is then not a problem that the properties of this proposed state do not match well with 
previously seen individual resonances. Nevertheless, it is rather consistent with an average of the lowest two states, 
similar to the p' case. 

Fourth, by adjusting the parameters of the a'^ states, the percent deviations of the Weinberg-type sum rules 
can change. The percent deviation of the second Weinberg-type sum rule is the most sensitive to such changes in 
parameters because the RHS of this sum rule is numerically very small. Each sum rule responds differently to each 
parameter. For example, the percent deviation of WSR-3 can be decreased by shifting the a'l mass to lower energies, 
however the percent deviation of all the other sum rules becomes worse. As an extreme case, the a'l can be shifted 
to lower energies accompanied by the reduction of the a'^ coupling so that WSR-2 and WSR-3 are both "precisely" 
satisfied; WSR-0 and WSR-1 are then violated by —2.1% and —3.6%, respectively. Although the overall deviations 
are smaller in this case than with the parameters in Table HI we believe that is it better to get the best possible 
agreement with the two sum rules which are most accurately know. 

Finally, the other fitted parameters seem to be in reasonable agreement with expectations. The ai mass is well 
within the PDG range [2^ , while its width is at the upper end of its range. The form- factor scales of around 1 GcV 
are of typical hadronic size, with the one for the ai somewhat, though not unreasonably, smaller. The threshold scale 
is consistent with decays into Att states for the p', while it is slightly larger than the physically expected Stt threshold 
for the a[. Although the coupling for the a[, ga[, seems large compared to the other states, its numerical value is 
dependent on the location of the continuum as discussed above. Overall, the chosen parameters result in spectral 
functions that fit the data and appear to be within reason of their physical interpretation. 

V. QCD SUM RULES 

Although the QCD sum rules were not used to initially constrain the model, it is interesting to determine to 
what extent the resulting spectral functions satisfy them. From Eqs. (|16p and (|17p. one can see that there are 
three condensates whose values are needed. The quark condensate will be set at the commonly used value, (qq) = 
(— 0.25GeV)'^. However, there is appreciable variability in the values of k and the gluon condensate in the literature. 
Their values will be tuned to optimize the agreement of the QCD sum rules given our input spectral functions. 

To quantify this agreement, we will use the method described in Refs. [itI, [T8|. In this method, one first defines 
a LHS and a RHS corresponding to the resonance part in the dispersion integral and the OPE with the continuum 
contribution from the dispersion integral subtracted, respectively. For our purposes, this division will be written as in 
Eqs. (dni) and PT)) with the continuum part of the spectral functions subtracted from both sides. Second, a range in the 
Borel mass, the so-called Borel window, needs to be established where one expects the best agreement between the two 
sides given their applicability limitations. We will define the lower limit, Mmin, in the standard way, i.e., as the Borel 
mass where the terms proportional to contribute at most 10% to the RHS, indicating that the OPE becomes 

unreliable at still smaller values of the Borel mass. For large Borel mass, the OPE is dominated by the term. 
However, this term is nearly trivially satisfied by the continuum. Therefore, in order to analyze the agreement due to 
the resonance states, an upper limit to the Borel window is introduced. This is typically chosen as the mass where 
the continuum contribution equals the one from the resonances. We have found that applying this procedure with our 
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FIG. 4: Graphical representation of the QCD sum rules as a function of Borel mass over the Borel window. The LHS (dispersion 
integral using the spectral functions of the previous section) of the sum rule is given by the solid curve while the RHS (OPE) 
is represented by the dashed curve. 

spectral functions results in a Borel window nearly double in size of what was previously found in the literature. This 
difference simply arises due to the significantly higher energies for the onset of our universal continuum. We therefore 
decided to redefine the upper limit of the window, Mmax, as the mass where the continuum contribution is half of 
the resonance contribution. This produces a Borel window of similar size as in previous studies. Instead of a 50/50 
split between continuum and resonances at the upper range, we thus have a 67/33 proportion between resonances 
and continuum. With these considerations, we found a Borel window of 0.85 GeV < M < 1.47 GeV for the vector 
channel and 0.89 GeV < A/ < 1.48 GeV for the axial-vector channel, which is comparable to earlier studies (T7l - [l9| . 

With spectral functions, condensate values, and the Borel window defined, the agreement between the two sides of 
the sum rule is measured by the value d defined as 



where AAf^ = Ml,„^ — M^-. One can think of d as the average deviation between the two sides over the Borel 

iritix mm o 

window. The optimization of the values for k, and the gluon condensate thus amounts to minimizing d. Optimizing 
the vector and axial- vector channels independently yields different values of k and the gluon condensate for the two 
channels. While the gluon condensate should be identical, the value of /t, representing correlations beyond the ground 
state, can, in principle, be different for different quantum numbers. For simplicity, we have decided to also assume 
a universal k value and to minimize the sum of the d values obtained from the two channels, that is, we minimize 
dtot = dv + dA- To estimate the uncertainty in the optimized values, the range in k and the gluon condensate which 
produces a combined deviation of less than 1% has been found. The minimization procedure gives k = 2.lj]Q'2 and 
i^Gl^) = 0.022 ± 0.002 GeV*. The optimized d values for vector and axial-vector channel are 0.24% and 0.56%, 
respectively. Increasing the upper limit of the Borel window, Mmax by 5% yields no significant change in the values of 
the parameters, but increases the deviation. Furthermore, when Ky and ka are decoupled, the optimization procedure 
yields rather similar results, (^GlJ = 0.022 GeV'', Ky = 2.1, ka = 2.0, dy = 0.24%, and dA = 0.55%, thereby 
justifying a posteriori the simplifying assumption of identical k values. 

Our values for k are well within the range previously found in the literature [H, |T3, [H, , while the value for 
the gluon condensate is in the upper range. For example, early applications in charmonium sum rules have extracted 
i^G^^) ~ 0.012 GeV^ 0- On the other hand, more recent studies of higher moments of the charmonium sum rules 
yield significantly larger values, of around 0.022 GeV^ [isl - liTj . close to our value. 

To illustrate the agreement of the QCD sum rules, we display in Fig. |4] both LHS and RHS for each channel, scaled 
by M^, as a function of M over the Borel window. With the sub-1% agreement, it is difhcult to distinguish the two 
curves in each figure. We note that an increase in k with fixed gluon condensate will shift the RHS of the vector 
channel down, while the RHS of the axial channel is shifted up (dashed curves), thereby making the agreement in 
both channels worse. Overall, the spectral functions constructed may be considered consistent with the QCD sum 
rules. 




(30) 
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VI. CONCLUSION 

In the present work, we have performed a combined analysis of Weinberg-type and QCD sum rules using vacuum 
spectral functions for vector and axial-vector channels constructed via quantitative fits to hadronic r-decay data. For 
the ground-state resonances, we employed a microscopic p spectral function and a Breit-Wigner ansatz for the ai. 
A novel feature of our approach is the introduction of an excited state in each channel which, in turn, allowed us to 
employ an universal perturbative continuum part at high energies. The universality of the continua is rather welcome 
in view of chiral degeneracy in perturbation theory, and leads to a higher onset energy than in previous works. While 
the excited vector state (p') was deduced from an accurate fit to data, we have found that the three lowest Weinberg- 
type sum rules can be quantitatively satisfied only if the existence of an a'l state is postulated. The latter lies outside 
the direct realm of the axial- vector r-decay data, but the extracted mass is compatible with (and in a sense confirms) 
an average of previously observed states. Furthermore, the resulting spectral functions were implemented into QCD 
sum rules, resulting in a sub-1% agreement with the operator product expansion (comparable to other state-of-the-art 
analyses). We believe that these spectral functions provide a good basis for future studies of medium modifications, 
to shed light on the long-standing problem of testing chiral symmetry restoration with dilepton data. 
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